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By a singular gauge transformation, the quasi-particle transport in the mixed state of high T c 
cuprates is mapped into charge-neutral composite Dirac fermion moving in short-range correlated 
random scalar and long-range correlated vector potential. A fully quantum mechanical approach 
to longitudinal and transverse thermal conductivities is presented. The semi-classical Volovik effect 
is presented in a quantum mechanical way. The quasi-particle scattering from the random magnetic 
field which was completely missed in all the previous semi-classical approaches is the dominant scat- 
tering mechanism at sufficient high magnetic field. The implications for experiments are discussed. 
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The general problem of quasi-particle transport in the 
mixed state of high T c cuprates is important, because si- 
multaneous measurements of thermal conductivities k xx 
and K xy provide a lot of information on the new physics 
pertinent to d wave superconductors. On the experimen- 
tal side, Krishana et al observed that in superconducting 
BSCCO and YBCO, at temperature T > 5K, the longi- 
tudinal thermal conductivity k xx {H) initially decreases 
with applied magnetic field H, then reaches a plateau 
(0. They also measured ^,|) thermal Hall conductivity 
K xy at T > 10K and extracted the thermal Hall angle 
taxi9 = ^ 2L . On the theoretical side, employing semi- 
classical approximation, Volovik pointed out that the cir- 
culating supercurrents around vortices induce Doppler 
energy shift to the quasi-particle spectrum, which leads 
to a finite density of states at the nodes . This effect 
( Volovik effect) has been employed to explain the above 
experimental observations of k xx by several authors ||. 
However, semi-classical method can not be used to calcu- 
late K xy . A fully quantum mechanical approach is needed 
to get any information on n xy . 

Starting from BCS Hamiltonian, Wang and MacDon- 
ald performed a first numerical calculation on quasi- 
particle spectrum in vortex lattice state By phe- 
nomenological scaling argument, Simon and Lee (SL) § 
proposed the approximate scaling forms for k xx and K xy 
for dirty d wave superconductors in the mixed state. An- 
derson H employed a singular gauge transformation to 
study quasi-particle dynamics in the mixed state. Franz 
and Tesanovic (FT) employed a different singular gauge 
transformation to map the quasi-particle in a square vor- 
tex lattice state to Dirac fermion moving in an effective 
periodic scalar and vector potential with zero average 
and studied the quasi-particle spectrum |J. Using FT 
transformation, Marinelli et al studied the spectrum in 
various kinds of vortex lattice states in great detail |Io| . 

In this paper, by considering carefully the gauge invari- 
ance overlooked by previous authors , we extend FT 
singular gauge transformation to include the curvature 
term which is important to K xy . After clarifying some 
important subtle points of the singular gauge transfor- 



mation, we prove exactly that there is no Landau level 
quantitation, then apply it to disordered vortex state with 
logarithmic interaction between vortices. We find that 
quasi-particle moving in the disordered vortex state of 
high T c cuprates can be mapped into charge-neutral com- 
posite Dirac fermion moving in short-range correlated 
random scalar potential and long-range correlated vector 
potential with zero average. The quasi-particle scattering 
from the long-range correlated internal gauge field dom- 
inates over those from the well-known Volovik effect and 
the non-magnetic scattering at sufficient high magnetic 
field H > H* ~ H c2 (^-) 2 . k xx satisfies scaling Eq. 13 



at H > H* , it approaches different pinning strength de- 
pendent plateaus due to the vertex correction shown in 
Fig.l. However, n xy satisfies scaling Eq.|l4], at H > H*, 
it increases as sfB.. 

We start from the d wave BCS Hamiltonian in the 
presence of external magnetic field: 



h + V(x) 



A 



At -h* - V(x) 



d(x) (1) 



with d-f (x) = Cf (x) , d[ (x) — c{ (x), V(x) is random chem- 
ical potential due to non-magnetic impurities. 

The non-magnetic impurity scattering part can be 
written as: 



H imp = V(x)d} a {x)T^d a {x) 



(2) 



We assume that V(x) satisfies Gaussian distribution 
with zero mean and variances A^ 

<V(x)>=0, < V{x)V{x') >= A A S 2 (x-x') (3) 

The gauge invariance dictates: 

h=±(p- e -Af-e F 
2m c 

A = [\Px ~ \d x <j),p y - -d y( j>}A(r) 
+ (Px - ^d x 4>)A(r)(py + ^d y <j)) 
+ (Py ~ o d y<t>) A (r)(Px + \d x <f) 



1 



A(r){p, 



M}] 



(4) 



where <f> is the phase of A(r) |l2] ]. 

The gauge invariance in A has not been taken into ac- 
count in Refs. |7|||. Although its correct treatment does 
not affect the linearized Hamiltonian ( see Eq.|] ), it is 
crucial to the curvature term ( see Eqj^ ) which is impor- 
tant to K xy . The similar gauge invariance was considered 
in Eq.2.12 of Rcf. |J). 

Following FT , we introduce composite fermion 
d c by performing a singular unitary transformation d = 
Ud r : 



H s = U-^HU, 



U = 



i4>A{r) 





i<pB (?) 



(5) 



where <j>A is the phase from the vortices in sublattice A 
and <f>B is the phase from the vortices in sublattice B. 

In this letter, we assume the thermal currents are suf- 
ficiently weak that the vortices remain pinned by non- 
magnetic impurities. Therefore, the transport properties 
of d is exactly the same as the composite fermions d c . It 
is easy to check that U commutes with r 3 , therefore the 
transformation leaves Eqj^ Hi mp invariant. 

Expanding H s around the node 1 where p — (p^, 0), we 
obtain H s = Hi + H c where the linearized Hamiltonian 
Hi is given by: 

Hi = v f (p x + a x )r 3 + v 2 (p v + a a )r 1 + v f v x (r) 
+ V(x)t 3 + (1 2,x^ y) (6) 



where v s 



\A is the total super- 



fluid momentum and a a = ^{v£ — v£) — h{V(f>A — V<^b) 
is the internal gauge field. Anderson's gauge choice is 
(j>A = 4>,4>b = or vice versa ||,[l5|]. 

We get the corresponding expression at node T and 2 
by changing Vf — > —Vf, v 2 — > — vi in the above Eq. 
The curvature term H c can be written as: 



Hr 



1 



{U a ,v a } 



n 2 



Ac 
2e, 



{n^njr 1 ] (7) 



Where II = p + a is the covariant derivative. H c takes 
the same form for all the four nodes. 

It is easy to see that v a (f) acts as a scalar scatter- 
ing potential, it respects time-reversal (T) symmetry, but 
breaks Particle-Hole (PH) symmetry |L7| . There are two 
very different kinds of internal gauge fields in Hf. V(x) 
is due to non-magnetic impurity scattering at zero field, 
a a is completely due to Aharonov and Bohm (AB) phase 
scattering ]l5[ from vortices generated by external mag- 
netic field. They both respect P-H symmetry. In general, 
V(x) breaks T symmetry, but T symmetry is restored 
in the unitary limit p8| . For general flux quantum a, 
a a breaks T symmetry, but T symmetry is restored at 
a = 1/2 Jl9| because a = 1/2 flux quantum is equivalent 



to a = —1/2 one due to the periodicity under a — > a+ 1. 
This is also the underlying physical reason why we are 
able to choose the two sublattices A and B freely with- 
out changing any physics. Due to this exact T symmetry, 
there is no Landau level quantization as claimed in S . 

In H c , the only term which breaks both P-H and T 
symmetry Q is ip^{p a , v a }ip = -iv a (V a i/' t V'-'0 t V aVO- 
This term will lead to thermal Hall conductivity to be 
discussed in the following p0| . 

From EqJ^, it is easy to identify the conserved charge 
currents at node 1: 



jlx = 1p{(x)VFT 3 l/jl(x), ji y = tl}\(x)V2T 1 tj}l[x) (8) 

with the currents at node 2 differing from the above ex- 
pressions by (1 — > 2, x — > y). It is known that the charge 
conductivity of d c corresponds to the spin conductivity 
of c electrons Because the spin o~ s and thermal conduc- 
tivities are related by Wiedemann- Franz law . For 
simplicity, in the following, we only give specific expres- 
sions for the spin conductivity which can be evaluated by 
Kubo formula. 

The Hamiltonian Hi + H c enjoys gauge symmetry 
U u (l) x U s (l), the first being uniform (or external ) and 
the second being staggered ( or internal ) gauge symme- 
try. Although the composite fermion d a is charge neutral 
to the external magnetic field, it carries charge 1 to the 
internal gauge field a a . 

We assume a randomly pinned vortex array with log- 
arithmic interaction between vortices. In the hydrody- 
namic limit, after averaging over all the possible positions 
of the vortices Ri, we find: 

< v a >=< a a >= 0, < v a (k)ap(k') >= 

< v a (k)vp(-k) >= n 2 (8 af3 - ^r-) k2 U ^ n 

<a a (k)a^-k)>=^(S a0 -^)^ (9) 

Where the vortex density is n v — U- — -ff— A . 

The first line in Eq.K is exact. The v — v and a — a cor- 
relators are the most general forms consistent with the 
incompressibility of the vortex system p2|| . The pinning 
strength will only enter as prefactors in front of n v . For 
notational simplicity, we suppress these prefactors. 

Because v and a are decoupled at quadratic order (see 
Eq.9), the long-range logarithmic interaction between 
vortices suppresses the fluctuation of superfluid veloc- 
ity, but does not affect the fluctuation of the internal 
gauge field. Therefore the scalar field acquires a " mass 
" determined by the vortex density, but the gauge field 
remains " massless " . The gauge field is a pure quan- 
tum mechanical effect, it was completely missed in all 
the previous semi-classical approaches Q,^] . Here we ex- 
plicitly demonstrate that being gapless, its fluctuation 
even dominates over the well-known Volovik effect in 
the low energy limit. It also dominates over that from 



2 



the non-magnetic scattering at sufficiently high field and 
low energy limit. 

In fact, in the weakly type II limit £ < A < d v , the 
superfluid velocity vanishes in the interior of supercon- 
ductor, the long-range correlated gauge potential in Eq.^ 
becomes the only scattering mechanism [ ^6| . 

From the Eqjp, it is easy to realize that the internode 
scattering ~ kp is weaker than the intra-node scattering 
~ Pq 2 by a factor a -1 = <SC 1. In the following we 
neglect the internode scattering. 

Up to the order of Gaussian cumulants, the scalar po- 
tential and vector potential are uncorrelated. However, 
they are correlated by Non-Gaussian cumulants. The 
lowest order non-Gaussian cumulants are the skewness: 

< v a (fi)vp(r2)a»y(ra) >=< a a (ri)ap(r 2 )a 1 (r 3 ) >= 

< Vaikijvp^^^kz) >= n 3 n v 5(ki + k 2 + k 3 ) 

-it a sku{k 2 ■ k 3 S f 3 7 - fcg-yfcg/a) 
(k\ + n v )(k\ + n v )(k\ + n v ) 

< v a {ki)ap(k 2 )a 1 (k 3 ) >= n 3 n v S(ki + k 2 + k 3 ) 
-ie aS k ls (k 2 ■ hS^ - k 2j k 3f3 ) 



(10) 



In fact, because any distribution function satisfies 
P[a a {x)] = P[—a a (x)], any correlators involving odd 
number of a a vanish. After coarse graining, the ex- 
act T symmetry of Hi is approximated by the average 
one. This approximation will lead to correct behaviors of 
self- averaging physical quantities. But it docs not apply 
to non-self-averaging quantity such as Hall conductance 
fluctuation. 

The discussion on k xx : The scalar potential capture 
the essential physics of Volovik effect: the quasi-particles 
energies are shifted by superfluid flow. Following the RG 
analysis in Ref. |l7| , it can be shown that the random 
scalar potential is marginally relevant, therefore gener- 
ates finite density of states at zero energy. Because 
kpl tr ~ kpd v S> ~ 5, the SCBA in standard im- 
purity scattering process can be applied to calculate the 
low energy scattering rate. In k, oj — > limit, it leads to 



VpTT n. 



d 2 p 



1 



(27T) 2 P 2 (p 2 +^)r2 



E 2 
p 



(11) 



where E% 



(vfki) + (v 2 k 2 ) , the retarded self-energy 
"E R (k, oj) = E(fc, iui — > u> + iS), the zero energy scattering 
rate is r = -7mS H (0,0). 

The above equation leads to the quasi-particle lifetime: 



iM~r ~Ao,/ — 



H 



(12) 



Now we look at the vertex correction to n xx due to the 
ladder diagram shown in Fig.l at T = 0. 



(^^) + (^T^) + = AA^^W 

1 

a b c 

Fig 1: Vertex correction to longitudinal conductivity 

Fig. la is just the bubble diagram. It leads to the well 



VfV 2 



known bubble conductivity |Tl],^l[]: a xx - 

By checking the integral equation satisfied by the ver- 
tex function T(p, iu>, iui + ifl) = r 3 (l + A(p, iu>, ito + iQ)) 
depicted in Fig.lc, we find A is at the order of 1 inde- 
pendent of n„, therefore o~ xx receives vertex correction 
of order 1 independent of n v |23] . This vertex correction 
was completely missed in the semiclassical treatments [B] . 

The random gauge field gives additional scattering 
mechanism, it has scaling dimension 2, therefore strongly 
relevant [H and dominates over Volovik effect at low 
energy limit. Similar SCBA to Eq.[ll] leads to logarith- 
mic divergent quasi-particle scattering rate 1/t/ which 
is not gauge-invariant anyway. But the vertex correc- 
tion similar to Fig.l removes the logarithmic divergence 
and leads to the finite gauge-invariant transport time 

The vertex correction to the bubble conductivity due 
to the non-magnetic impurity scattering Eq.[s] among the 
four nodes was calculated in Ref. |ll]], it was found to 
be negligible. It is obvious that the two vertex correc- 
tions are different due to different scattering mechanisms, 
therefore the two conductivity values are different, al- 
though the bubble results are the same. 

At finite temperature, the T dependence comes solely 
from the Fermi function, a xx should satisfy the following 
scaling ( T c ~ A ): 



a xx (H,T) = F^aTrtr) = F x (a^^) 



(13) 



This scaling is consistent with Simon and Lee [Q us- 
ing phenomenological scaling argument. Our derivation 
bring out explicitly the underlying physical process: the 
quasiparticle scattering due to the long-range correlated 
random gauge potential. Pushing further, we conclude 
that in the high field limit H ^> a 2 H c2 (^-) 2 , a xx should 
approach the T = value F\ (0) = a xx (0) at the order of 
1. This value depends on not only the anisotropy param- 
eter a = Vf/i>2, but also the pinning strength appear- 
ing in Eq.^. This dependence could explain the different 
plateau values observed in the experiments pj . 

Unfortunately, the value a xx (0) is hard to be sorted out 
experimently due to the large background contributions 
from phonons [jlpl . 

We start with Hi. In order 



The discussion on k 



to get a non- vanishing n xy , we must identify terms which 
break both T and P-H symmetry. As shown previously, 
Hi respects exact T symmetry, therefore < a xy >= to 
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the linear order. We have to go to the curvature term 
Eq.^ to see its contribution to < a xy > . 

The first contribution comes from the skewness 
between Volovik term and the {p a ,v a } term < 
dv a (fi)v/3(f2)vj(f2) > in Eq.flO| which breaks both T and 
P-H symmetry. 

Just like in k xx , the skewness between scalar and ran- 
dom gauge field < dv a (fi)af3(r2)a 7 (f 3 ) > in Eq.|lO| gives 
additional scattering mechanism. It even dominates over 
the pure scalar skewness at low energy limit. 

Due to the antisymmetric tensor e a s, we find o~2 X y — 
—o~i yx = &i X y. Because both skewnesses are even under 
Vf — > —Vf,V2 — ► V2, it is easy to find that o\ xy — o\ xy , 
therefore a xy = a lxy + a 2xy + o\ X y + a^xy = A&ixy 

Because the {p a ,v a } term contains one more deriva- 
tive, simple power counting leads to 



This scaling is consistent with Simon and Lee |?]] 
by phenomenological scaling argument. Our derivation 
bring out explicitly the leading contributions from the 
{Pa,v a } term in the curvature term and also its small 
numerical factor Pushing further, we conclude that 

in high field limit H ^> b 2 H c2 (^-) 2 , cr xy should in- 
crease with H as ^yT^F 2 (0). Taking a ~ 10, H ~ 

10T, H c2 ~ 150T, we find the prefactor is about 1/40, so 
K xy is smaller than k xx by a factor of 1/40. The smallness 
of K xy make it difficult to measure experimentally. 

The most recent data at 10K < T < 30K for n xy /T 2 
H was shown to satisfy quite well the scaling n xy /T 2 = 
F(by/H /T) which follows from Eq|l4|, however, it con- 
tinues to decrease up to H = 14T instead of increasing 
linearly with y/H/T as follows from Eq.|l4]. The discrep- 
ancy may be due to the inelastic scattering at T > 10K 
not considered in this paper. For technical reason, so far 
the data is not available at T < 1QK. 

In conclusion, we point out that the dominant scatter- 
ing mechanism at sufficient high magnetic field is due to 
the long-range correlated random gauge potential instead 
of the well-known Volovik effect. 
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